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TRANS I ENT- FORCED CONVECTION FILM BOILING 
ON AN ISOTHERMAL FLAT PLATE 


INTRODUCTION 


Under the starting conditions of the pumps to feed cryogenic fluids 
for combustion in the thrust chambers of space vehicles, the temperature 
associated with fluid passages causes excessive vapor formation caused by the 
large temperature differences between the cryogenic fluid and the surrounding 
walls. This sudden vapor formation often reflects in the blockage of flow 
referred to as "vapor lock. " In turn, the NPSH of the pump will be drasti- 
cally affected. For remedies to the problem, the following methods are 
considered: 

1. Preconditioning the flow passages [1] before starting the pump. 

This method has been used in the missions undertaken so far. 

2. The possibility of using a pump capable of handling two-phase 
fluids revealed that a suitable design [2] of the blades can handle vapor liquid 
mixture ratios up to 40 percent without seriously impairing the capabilities of 
the pump. 

3. The circulation of the cryogenic fluid with increased velocity until 
the walls attain the desired temperature for a normal operation of pumps. 

For space vehicles operating under low-gravity conditions for long 
durations and necessitating frequent restarts, methods (2) and (3) are being 
considered. It is evident that two-phase pumping cannot be used especially 
when the volume of vapor generated is large. Thus, forced convection is 
considered as a promising remedy. 

The high temperature difference between the wall passages and the 
cryogenic fluids results in film boiling during the initial phase of "chill down. " 

The available literature on forced convection film boiling [3, 4, 5, 6] is 
centered around the steady state analysis while no investigation of the transient 
case is known to the author. However, the method of characteristics has 
been used to analyze the transient flow and heat transfer for forced convection 
in the entrance region of flat ducts [7] and for free convection on a vertical flat 
plate [8] . In the present work, transient -forced convection film boiling on 



a flat plate is analyzed by similarity and local similarity transformation tech- 
niques. For making the problem tractable, the isothermal condition of the 
plate is assumed. When the external heat which soaks into the system is 
balanced by the heat carried away by the flowing vapor-liquid system, this 
assumption may be justifiable . 

PHYSICAL MODEL 

Figure i shows the physical model and the coordinate system used in 
the analysis. The plate is assumed to be at a uniform temperature, T , 

above the saturation temperature of the fluid. The liquid flows over the plate 
with a constant velocity, , before touching the plate and the velocity, 

U , far away from the plate remains unaltered throughout. The liquid in the 

OO 

immediate neighborhood of the heated plate vaporizes and a thin layer of 
vapor is formed. The formation and growth of this layer is assumed to be 
continuous and stable. The other assumptions include the following: 

1. The liquid is at a saturation temperature far away from the plate. 

2. Effects of radiation and viscous dissipation are neglected. 

3. Two-dimensional laminar conditions prevail. 

4. All fluid properties are invariant. 

5. The order of magnitude of velocity and temperatures in the 
direction of flow merit consideration of boundary layer approximations. 


FORMULATION OF THE PROBLEM 


Under the above set of assumptions, the equations of mass, momentum, 
and energy are: 

1. Vapor Layer 


8X + 8Y 


!IL + U au 8u efu 

8t t 8X 8Y 8Y 2 


2 



and 


+ u — 

at* ax 


+ v 


9T 

a y 
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a 2 T 

8Y 2 


and 


2. Liquid Layer 


au au 

— — + — - = o 

ax aY 


9U i * d A v Oh 

at* + U i ax + V i aY v i aY 2 


Consideration of the energy equation for liquid layer is not necessary since the 
liquid is assumed to be at a saturation temperature. The following nondimen- 
sionalizing variables are used: 


and 


U 



OO 


V 



OO 


X = 


X 
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Re 


U L 

OO 


V 



Pe = Re Pr , 



da) 
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where 6 is the boundary layer thickness. Similar transformation variables 
for liquid layer are: 



OO 



U 

OO 


and 


U L 

OO 



(lb) 


The resulting equations are: 

1. Vapor Layer 

8u 8v . 

— + — = 0 , 

9x 9y 


(2a) 


8u 8u du. _1_ a 2 u 

at + u ax + v 3y Re ay 2 

and 

aT _ay aT l a 2 T 

at ax ay Pe ay 5 


(2b) 


(2c) 
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2. Liquid Layer 
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The boundary and initial conditions for the problem are: 

tsO u = v = 0 everywhere 

T = T at y = 0 
w 

T = T for y > 0 

oo J 

t>0 u=v=0;T=T aty=0 

w 

T = T at y = A 

CO 

u„ 1 as y — 

f J 

where A is the vapor layer thickness. 


(2d) 


(2e) 


(3) 


ANALYSIS 


For impulsive start of motion on a flat plate, when there is no heat 
transfer the problem can be solved as an exact solution of Navier Stokes 
equations [9] by introducing similarity transformations. The growth of 
boundary layer and the flow are independent of the axial coordinate. 

For forced convection film boiling under steady state conditions, 
similarity transformations have been found [3] in the two-phase boundary 
layer analysis. The solutions are independent of time and are functions of x. 
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In an attempt to develop the similarity transformations to describe the 
complete transient process the following variables are sought. The stream 
functions and are defined as: 

i . Vapor Layer 


and 



v = - 


9 * 

9x 


(4a) 


2. Liquid Layer 
9 


u f 9y 


and 


__ 

V i 9x 


(4b) 


satisfy the continuity equations (2a) and (2d) . The next set of variables used 
are 




A(x,t) ’ 


*(x,y,t) = A (x, t) f(T 7 ) , 


and 

T(x, y, t) - Too = 
(T w -Too ) 

and for the liquid layer, 


H(r)) , 


(5a) 
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__z 

A £ (x,t) 


and 


(5b) 


(x,y,t) = A £ ( x, t) F(r?) , 


where 1 


and 


A(x,t) = \[2 Re ^ (Cjx + C 2 t) 1 / 2 


A £ (x,t) = 42 Re £ ^ (Cjx + C 2 t) 


(5c) 


( 5d) 


Introducing these variables in equations ( 2) , the following set of ordinary 
differential equations are obtained: 


f'" + C 2 rj f" + Ciff” - 0 , 

H" + Pr (C 2 Ht)+ Ctf H') = 0 


(6a) 

(6b) 


and 


F ,m + C 2 tj F" + CiFF" = 0 . (6c) 

It should be noted that Cj and C 2 are considered initially as constants 
with respect to the independent variables x, y, and t. Referring to equation 
(5c) the variation of these parameters could be estimated. The value of zero 
for Cj indicates that the solutions are independent of x and represent the 
initial phase of the transient growth, further, the value of C 2 for Ci = 0 can 
be estimated a priori from the available solutions [9] as about 8. 0. However, 
this value will be obtained from the present solutions also. For C 2 = 0, the 
solutions are independent of time. Hence, for steady state, C 2 = 0. Thus, 
the range of C 2 will be from 0 to about 8.0. Similarly, the range of values of 


i. The method of obtaining these transformations are given in Appendix A. 
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Cj can be estimated. For a given value of AT/h^Pr , there is a 

unique value of under steady state conditions. Thus, variation of will 

be from 0 to b where b depends on the parameter C AT/h„ Pr. The 

P fg 

boundary and initial conditions in terms of transformed variables become, 

77 = 0, f = f* = 0 H = i 
77 = 1, H = 0 

Vg *■ 00 > F'-*l (7) 

INTERFACE MATCHING CONDITIONS 

At the vapor-liquid interface, the velocity, shear stress, mass, and 
energy balance conditions should be satisfied. An analysis similar to steady 
state case [3] is carried out. The result of this interface matching under 
boundary layer assumptions is as follows: 

1. Vapor liquid interface is represented by 77 = 1 or 77 ^ = R t , where 

Rl= (v/^) l/2 . (8a) 

2. The velocity balance yields, 

F r (Ri) = f' (1) . (8b) 

3. Shear stress matching gives 

F” (Rj) = Rf" (1) , 

where 

R = [W)/(PM)J 1/2 . (8c) 
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4. The continuity requirements demand, 


F(R 1 ) = Rf (i) 


(8d) 


5. The energy balance can be written at the interface as, 



The left-hand side denotes the heat transferred by conduction at the interface 
and the right-hand side represents the heat convected during vaporization of 
the liquid at the interface. The subcooling of the liquid is neglected. In 
terms of the transformed variables, equation (8e) becomes, 


AT 

P _ C 1 f (1) 

h, Pr -H’ (1) 
fg 


Thus, a unique relation exists between the constants, and C 2 , and the 

property function, (C AT/h. Pr) . It can be noted that C x plays the same 

P 

role as [3] in steady state conditions. 

Since equations ( 6) contain rj, the value of r?^ at the interface cannot 

be taken as zero as in the steady state case [3] . Hence, the relation 
[equation (8a)] should be used. Hence, a new parameter, R lf also 
appears apart from other parameters; i.e. , Pr and R. 


Steady State Heat Transfer 

Under the steady state conditions, the rate of heat transfer becomes, 


q(x) = - K 


/ 9T(X,Y) 
\ 3 Y 


■) 


Y=0 


(9a) 


In terms of the transformed variables using C 2 = 0 for steady state, 


9 


q(x) = Re 1 / 2 

L \J 2C t x 


(9b) 


Defining the local heat transfer coefficient and Nusselt number as 


h = sM. 

x AT 


and 


h X 

XT X 

Nu = — — — 
x K 


(9c) 


the expression for heat transfer relation can be: 


Nu 


Re i / 2 
x 


(y- 


-R 


\[2C[ 


H s (°) 


( 9d) 


Fluid Flow Relation 


The shear stress at the plate is 


t = M 


9U\ 

3Y 'y=o ' 


(10a) 


In terms of the present variables under steady state, equation (10a) can be: 



He */ 2 = f"(0) 

X n/2C7 


(10b) 
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Transient Parameters 

The rate of heat transfer from the plate is: 


q(X,t 1 ) - -K 


BTtX.Y.tj) 

8Y 



(lia) 


In terms of the present variables, equation ( ila) becomes: 


Q(x,t). = 


-KAT H tr (0) 


(lib) 


The Nusselt number in the transient region is: 


Mn — q(x > ^ ^ — - h* ( oi x 
Nu tr- ATK ' H tr ° A 


(11c) 


Under the steady state conditions, the Nusselt number is given in equation 

a s 


( 9d) . The ratio of Nu^ r to Nu R will be useful for comparison purposes: 


Nu, 


tr 


H’ (0) 
tr 


Nu s H s (0) U 


(lid) 


A similar expression for shear stress can also be obtained. 


Total Heat Transferred During the Transient 

The total heat transferred until steady state is: 
t S 

q tr (x) = / q(x,t) dt , (12a) 


where tg is the time necessary to reach steady state. Equation (12a) , in 
terms of the present variables, becomes: 
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t S H' (0) 

-KAT r tr ( 1 „ 

V (X) = I” J A dt 


(12b) 


Defining an average rate of heat transfer during transient, equation ( 12b) 
becomes: 


- . . 1 / KAT\ r 

v (x) = -r \—)l 


*s h’ (0) 

- tr 


dt 


(12c) 


The ratio of (x) to the rate of heat transfer at steady state is important 

in presenting the main features of transient heat transfer. The ratio, after 
simplification using equation (9b) , becomes: 


V< x) = r 1 H 'tr (0) 
Q(X) 0 A H g ' (0) 


(12d) 


Evaporation Rate Expressions 

The rate of vapor carried away by flowing fluid is: 
5 

M (X,t t ) = f p UdY 
C 0 


M ( x, t) = pU L f ud; 


(13a) 


- (pU L)A f (1) 
°° tr 


The total vapor, carried away during transient is: 


M (x) = (pU L) / A f ( 1) dt . 


co J tr 


I 


Rate of evaporation at steady state is: 


(PU^L) 

M (x) = — n/ 2C iX f c (i) . (13b) 

Re 1 / 2 S 


The average rate of vapor convected during transient is: 

, *B 

M (x) = (pU L) r / 41 (1) dt . 

c l s 0 1 


(13c) 


The ratio of the average transient convection of vapor to the evaporation rate 
at steady state is: 


M (x) 
c 

M(x) 


i 


/ 


f s (i) 



(13d) 


Apart from the vapor convected away by the fluid, there is accumula- 
tion of vapor also which increases the boundary layer thickness. The total 
accumulation of vapor at any distance, x, until steady state is: 

X 6 

M (X) = p / dX / dY 
a 0 0 


M (x) = p L 2 Re 
a 


- 1/2 


x 

f \l 2C lX 

0 


dx 


( 13e) 


= 0.66 \J 2C lX Re"" 1/72 (pL 2 ) x 

The average rate of accumulation is: 

_ M (x) _ l/2 

M (x) = — = 0.66 \l 2CjX Re / pU L 

a 



( 13f) 
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The ratio of average rate of accumulation during transient to the steady state 
evaporation rate is: 


M (x) 
a 

M(x) 


0.66 

fgd) (t s /x) 


(i3g) 


The evaporation of vapor will be the sum of vapor accumulated and vapor 
convected. The ratio of the average rate of evaporation during transient to 
the steady state value will be the summation, 

M. (x) M (x) 

tr _ e 

M(x) “ M(x) 


M (x) 
a 

M(x) 


(13h) 


SOLUTIONS 


The set of equations (6) with boundary conditions (7) and interface 
matching relations ( 8) has been solved numerically on a digital computer for 
parametric values of Cj [0(0. 2)1. 5], C 2 [0( 0. 5) 2( 1) 8] , R(0.1,0.01, 0.001), 
R t ( 10, 25, 100) , andPr(0.5, 0.65, 1.0). These ranges were selected through 
an a priori knowledge from steady state results [4] to yield results of 
practical importance. The detail of the method of numerical solution is the 
same as that used for steady state analysis [6] . 


DISCUSSION OF RESULTS 


The numerical solutions of the set of equations (6) yield different 
values of R^ for each set of c i and C 2 , Pr > and R i- Representative 

results for R=0.1, Pr=i, and R 1 =25.0 are illustrated in Figure 2 as a 
plot of R^ versus Cj for parametric values of C 2 . It can be noted that the 

curve C 2=0 represents the unique relation between R^ and Cj [equation 

(8f)] under steady state conditions. This curve is the same as the one 
obtained by Ito and Nishikawa [4] for boiling without subcooling. The curves 
for finite values of C 2 indicate the transient period culminating in the starting 
phases of motion for C^O, the ordinate. 
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Initial Phases of Motion 


The above results are replotted in Figure 3 with Cj and C 2 on the 
axes and R^ as the parameter. It can be observed that for small values of 

and C 2 the curves will be straight lines. Also, small values of C 2 
correspond to the condition near steady state while small values of Cj 
indicate the initial phases. By using the equation for a straight line in the 
initial phases as 


C 2 - -a* Cj + a 2 


(14a) 


where a t is the slope of the straight line and a 2 is the intercept on the 
ordinate, the expression for A can be put as: 


A 

\J 2a 2 t 



(14b) 


For a value of — — = a t , the boundary layer thickness becomes equal to 1 
t 

irrespective of the value of and the variation of A depends on only 

a 2 

time. Thus, the value of a t indicates when the effect of x comes into the 

x 

picture. Further, for all values of - — > 1, the one -dimensional assumption 

L Si J 

will be made and the solution will be : 


A = \l 2a 2 t . ( 15a) 


It should be noted that for x/t > a 1# the value of A A J 2a 2 t becomes > 1, 
resulting in an unrealistic situation where the boundary layer thickness 
becomes greater than that for x/t = aj. 


Motion Near Steady State 

A similar analysis for small values of C 2 can now be made using 
Cj = -b 1 C 2 + b 2 as the equation for the straight line. 
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A 


A 


(15b) 


\ / 2b 2 x 



As before, the value of t/x = bj indicates when the effect of t can be 
neglected. For values of (t/b 1 x)>i, the results give rise to steady state 
solutions. It can be noted, as before, that for values of t/x > b 1} the value 
of A becomes >i and increases monotonically. Hence, the assumption that 
A = 1 for all (t/bjx) > 1 will be made. 

Thus, the slopes of the straight lines at the two ends of the curves in 

Figure 3 give the demarcating values for t/x. Figure 4 brings about these 

results clearly. The plot of t versus x consists of two straight lines having 

slopes i/aj and b 1# respectively. For all values of t below the line OA, 

the solution will be given by the one -dimensional heat conduction equation. 

Similarly, for the region above the line OB, the steady state solutions apply. 

The region in between these two is called the intermediate region. Similar 

classification was observed in the transient analysis of free convection [8] . 

However, in forced convection [7] , only two regimes were distinguished 

since fully developed velocity profiles were used in the analysis. The slopes 

a, and b< depend on the parameters R, R, , and Pr; the effect of R< 

h 

being negligible. The steady state time (t„) and the time (t ), below 

o o 

which one -dimensional conduction is predominant, depend on the parameters 

R, R, , and Pr which is shown in Figure 5. The ratio (t/x) is plotted 
h 

against R^ for R = 0. 1 and 0.001 and for a Pr = 1.0. The steady state 

value (t^/x) increases with increase in the value of R, up to about 1. 5 
o h 

and then becomes almost horizontal. This indicates that for vapors having 

large heat content (C AT) compared to the latent heat of vaporization, the 

P 

steady state time becomes almost constant. However, the value of (t Q /x) 

increases slowly with R^. The effect of R is shown by plotting the curves 

for R = 0. 001. The effect is less than 1 percent for all values of R^. Thus, 

the interfacial shear has negligible effect on transient film boiling in deciding 
the steady state time. A general observation of the figure is that the 
intermediate region occupies about 85 percent of the entire transient 
development. 
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Hitherto, all efforts were directed towards analyzing the problem 
[7, 8] for very small intervals of time and the study near the steady state. 

Such analyses did give the time to achieve steady state and the time indicating 
when the flow can be considered one dimensional. But, the method of 
characteristics used therein did not bring about the features of the intermediate 
regime . 


Intermediate Regime and the Local Similarity 


It is seen already that the ratio, C 1 /C 2 , represents a transient 

parameter. To obtain the development of boundary Jayer thickness in the 

complete transient regime, Figure 6 is plotted with A as the ordinate and 

(t/x) as the abscissa. The construction for the boundary layer growth is 

illustrated for R. = 1.0, Pr = 1.0, and R = 0.1. Using Figure 3, the 
h 

values of Cj and C 2 are read off for R^ = 1. Knowing the value of Cj 
and C 2 , A can be plotted with respect to (t/x) by the expression: 



(16) 


where b 2 is the steady state value of (= 1. 56 for R = 1.0). For 

several combinations of values of and C 2 , a family of curves 1.2. ... 7 
are plotted. It is obvious that none of the single curves satisfy the boundary 
conditions. This means that there is no similarity transformation applicable 
to the intermediate region. Figures 7, 8, and 9 indicate that the envelope of 
these curves satisfies the boundary conditions. This means that for each 
(t/x) , the values of C* and C 2 have to be changed to yield the necessary 
solution. Thus, by constructing similarity transformations that are locally 
valid, the solutions can be sought. 

It was observed during the process of obtaining A that at a fixed 

value of (t/x) , a minimum of A exists and this value A . is in fact the 

min 

enveloping curve obtained in Figures 7, 8, and 9. For obtaining this 

minimum value, equation ( 13) is used. For A . , 

min 


8 A 
8C t 


= 0 
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which gives 


(t/x) 



( 17 ) 


Thus, a relation exists between the slope of the plots in Figure 3 with (t/x) . 
Similar expressions giving the values of (t /x) and (t c /x) were already 

O o 

discussed. Now, the process of obtaining the solution can be summarized as 
follows: At any given time (the distance and properties of the fluid) the slope 

^ is found from equation (14) . The particular combination of Cj and 

C 2 for a given fluid can be obtained from Figure 3 and the value of A is 
calculated from equation ( 13) . By repeating the process for different (t/x) , 
the entire boundary layer growth can be obtained. Thus, the principle of the 
local similarity is used to throw more light on the aspects of intermediate 
regime . 

The concept of local similarity is used extensively in the boundary 
layer theory [10] . A recent discussion of this aspect can be found in Reference 
11. The essential features of this concept are that: 

1. The changes on the time scale are sufficiently small so that the 
time derivatives in the transformed governing equations can be neglected. 

2. The time variable, t, can be treated as a parameter wherever it 
occurs explicitly. 

In the present case, for the intermediate regime, Cj and C 2 become 
functions of t and, hence, they can be treated as parameters in solving the 
governing equations. 

As mentioned in Reference 11, no criteria for establishing the validity 
of the local similarity method appear, except that the time derivative in the 
transformed differential equations be small. Again, the smallness is only 
qualitative and no a priori decision has been formulated. Thus, under the 
assumptions of local similarity, the entire boundary layer growth has been 
obtained (Fig. 6) for a representative value of R, R^, and Pr. The 

"patching together" of local solutions [10], a characteristic feature of local 
similarity is clearly demonstrated while constructing the curves representing 
the boundary layer growth in Figure 6. Similar curves for different values of 
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R^ and Pr can be obtained. The variation of A for different values of R^ 
is shown plotted in Figure 10 using a nondimensional ratio of times ( t/t_) as 

O 

the abscissa. This mode of using the ratio of transient time to steady state 

time as the abscissa suggests that only slight variations (maximum of 5 

percent) exists between the several curves obtained for different R, , and a 

h 

mean curve could be plotted as shown. 

Figures 7, 8 , 9, 11, and 12 present the transient developments of 
velocity and temperature profiles. In Figure 11, f" is plotted against the 
independent variable, T), for the vapor region and (17 - Rj) for the liquid 

region. It was seen in Figure 3 that there are several combinations of C A 
and C 2 which give rise to the same value of R^. Since the value Cj = 0 

represents the case of pure forced convection, the ratio C^/C 2 was chosen 

for presenting the results. Thus, C x /C 2 = °° represents the steady state 

solution and the positive values denote the transient conditions. Figure 11 

shows one of the representative sets of velocity profiles for R =1.0 and 

h 

Pr = 1.0. For Ci/C 2 = 0, the growth of velocity profile is the steepest and 
coincides with the results of Reference 9. Similarly, the steady state results, 
Cj/Cz = °°, are identical with those of References 4 and 5. The transient 
development of velocity profiles is clearly shown in the gradually flattening 
velocity profiles for C l /C 2 values = 0. 015 to 0. 82. The effect of interfacial 
shear is brought about by the plot for R = 0. 001. Except very near the 
steady state, the effect appears negligible. The effect of R 1 is less than 
0. 1 percent in all cases. The corresponding transient developments of 
temperature profiles are shown in Figure 12. 

Figure 7 shows the effect of R^ on the velocity profiles near steady 

state and near (t Q /x) represented by values of = 1. 0 and 0. 2. One of the 

interesting observations is that all the profiles asymptotically reach the forced 
convection profile. For small values of R^, the velocity profile becomes 

almost linear in the vapor and liquid regimes. Similar features for tempera- 
ture profiles are shown in Figure 8 . 

The effect of the Prandtl number on the temperature profiles is 
demonstrated in Figure 9. Only the extreme cases are illustrated for clarity. 
For C 1 /C 2 = 0 (forced convection) , the effect of decreasing the Prandtl 
number is to decrease the boundary layer thickness. For C 1 /C 2 = , the 
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the decrease in Pr makes the profile linear. For values of C^/Ci in 
between these extremes, similar observations can be made. To illustrate the 
results of the local heat transfer parameter and the rate of evaporation during 
transient, the variation of H 1 (0) and f( 1) is studied for the parametric 
values of Cj C 2 , R, and Pr. Representative results are shown in Figures 
13 and 14. The variation of H* (0) with Cj for different values of C 2 is 
illustrated in Figure 13. The plot clearly shows a linear variation of H 1 (0) 
with Cj for each C 2 . With this family of straight lines and by using 
Figure 3, the variation of H 1 (0) with Cj for various values of R^ can be 

obtained. Such curves for representative values of R^ = 0.1, 1.0, and 2.5 

are shown which indicate a decreasing trend with respect to Cj, similar to 
the variation of C 2 in Figure 3. A study of these curves for different values 
of R shows that the effect is less than 1 percent in all cases except for steady 
state case (C 2 =0) . However, similar curves that result for various values 
of Pr are not shown here . 

Figure 14 shows the same characteristic features as Figure 13. 
Further, the effect of interfacial shear on f( 1) is also less than 1 percent 
except for C 2 =0. These curves vary for different Prandtl numbers because 
R^ depends on Pr. 

The ratio of the local heat transfer parameters at any instant of time 
to that at steady state [equation ( lid)] , as a function of t/t c , is plotted in 

Figure 15 for different values of R^. For t= 0, the heat transfer coefficient 

goes to °o in view of the sudden temperature increase of the plate. Through 

the transient process, the heat transfer coefficient remains higher than the 

steady. state values as seen in Figure 15. For fluids of lower R, values, the 

h 

heat transfer coefficient decreases more rapidly than for the fluids with 
higher values of R^. 

The average rates of heat transfer and evaporation during transient 

as a function of R, are shown in Figure lb. The plots of equations ( 12d) 
h 

and ( I3e) indicate that the average heat transfer rate during the transient 
varies from 1. 5 to 2 times the steady state heat transfer rates for R^ 

values of 0. 1 to 2. 5. Thus, the heat transfer rate in transient is much higher 

compared to that during steady state. Further, the rate of heat transfer 

increases for increasing values of R . The plots of average rates of total 

h 

evaporation, vapor convected away, and vapor accumulated presents the 
following: 
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1. The vapor convected away is 10 to 15 percent lower than that at 
steady state. 

2. Most of the vapor formed is consumed by the developing boundary 

layer . 


3. The average rate of vapor convected and accumulated decrease 
with increasing values of R . The rate of decrease of M /M is very large 

n. 3. St 

compared to M /M . 

C St 

4. The average rate of evaporation during transient is 2 to 15 times 
that at steady state. 


CONCLUSIONS 


A new approach using similarity and local similarity concepts has been 
adopted to obtain solutions for the boundary layer equations for transient -forced 
convection film boiling on an isothermal flat plate. 


The entire transient regime could be classified into three regions: the 
one -dimensional conduction regime, intermediate region, and steady state 
region. Exact solutions are obtained for the initial and final phases of motion 
and the expressions for times demarcating these regions. Using the local 
similarity concept, the intermediate phase is analyzed. 


The effect of interfacial shear is less than 1 percent in determining the 
steady state time. The parameter (C AT/h Pr) plays a very important 

role in deciding the rate of heat transfer and evaporation. Larger values of 

(C AT/h„ Pr) increase heat transfer and decrease evaporation. 

P fg 


The complete history of boundary layer growth (until steady state 
conditions) is obtained from the local similarity analysis. When the ratio of 
time at any instant to the steady state time is used, the boundary layer 
profiles for several values of R collapse to within 5 percent of a mean curve. 


The average rate of heat transfer during transient is 50 to 100 percent 
greater than the steady state heat transfer rates while the average evaporation 
rate is 2 to 15 times larger. The rate of accumulation of vapor to form the 
boundary layer gives the major cpntribution to the total evaporation. It varies 
from 1 to 14 times the steady state value. The rate of vapor convected away 
is about 10 to 15 percent smaller than the steady state values. 
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Figure 4. Classification of different regimes 
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Figure 16. Average heat transferred and evaporation durin 









APPEND IX A 

INTRODUCTION OF TRANSFORMATIONS 


When the transformations in equations (5a) and (5b) are introduced into 
the original set of equations ( 2) , the following partial differential equations 
result: 1 


f 

777717 


+ ARe (A T) f 

t VV 


+ A f f 

X 7)7) 


= 0 


H + A Pe ( A, 7) H + A f H ) = 0 , 

7)7) t 7) X 7)' 


and 


f + A. Re. (A. . 7) F +A . FF ) = 0 

7)7)7) f i t\l 7)7) X,i 7)7)' 


It can be seen that the above equations become ordinary ones if AA^ 

and A A can be made constants. These conditions resulting in two 
x 

simultaneous partial differential equations can be solved to yield. 


A = (constant) (CiX+ C 2 t) 


i/2 


where again C* and C 2 are arbitrary constants. The constant outside the 

-1/2 

parenthesis is chosen as n 2 Re so that Re does not appear explicitly 
in the ordinary equations. For the liquid layer, 


A l = n/T Re^ 1/2 (C lX + Cgt) 1 / 2 


i. Subscripts 7), x, and t are used to denote partial differentiation with 

respect to those variables. 
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